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Abstract 

The effective bottom Yukawa couplings are analyzed for the minimal super symmetric ex- 
tension of the Standard Model at two-loop accuracy within SUSY-QCD. They include the 
resummation of the dominant corrections for large values of tg/3. In particular the two- loop 
SUSY-QCD corrections to the leading SUSY-QCD and top-induced SUSY-electroweak 
contributions are addressed. The residual theoretical uncertainties range at the per-cent 
level. 



1 Introduction 

The Standard Model (SM) predicts the existence of one scalar Higgs boson which constitutes 
the remainder of electroweak symmetry breaking by means of the Higgs mechanism [ij. In all 
experiments this particle has escaped detection so far. Due to the hierarchy problem in the 
context of Grand Unified Theories supersymmetric (SUSY) extensions of the SM are considered 
as the most attractive solutions (if. The minimal supersymmetric extension of the SM (MSSM) 
requires the existence of five elementary Higgs bosons, two neutral CP-even (scalar) bosons h, 
H, one neutral CP-odd (pseudoscalar) boson A and two charged bosons if^. At lowest order all 
couplings and masses of the MSSM Higgs sector are fixed by two independent input parameters, 
which are generally chosen as tg/3 = V2/V1, the ratio of the two vacuum expectation values 
Vi^2, and the pseudoscalar Higgs mass Ma- Including the one-loop and dominant two-loop 
corrections the upper bound on the light scalar Higgs mass is Mh < 135 GeV [3|. More recent 
first three-loop results confirm this upper bound within less than 1 GeV The couplings of 
the various Higgs bosons to fermions and gauge bosons depend on mixing angles a and /3, which 
are defined by diagonalizing the neutral and charged Higgs mass matrices. They are collected 
in Table [1] relative to the SM Higgs couplings. For large values of tg/3 the down-type Yukawa 
couplings are strongly enhanced, while the up-type Yukawa couplings are suppressed, unless the 
light (heavy) scalar Higgs mass ranges at its upper (lower) bound, where the couplings become 
SM-like. This feature causes the dominance of bottom- Yukawa-coupling induced processes for 
large values of tg/3 at present and future colliders as Higgs decays into bottom quarks and 
Higgs bremsstrahlung off bottom quarks at hadron and e~^e~ colliders. Moreover, Higgs boson 
production via gluon fusion gg — h,H,A is dominated by the bottom-loop contributions for 
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Table 1: MSSM Higgs couplings to SM particles relative to SM Higgs couplings. 



large tg/3. Thus, the bottom Yukawa coupling determines the profile of the MSSM Higgs bosons 
for large tg/3 to a large extent. 

The negative direct searches for the MSSM Higgs bosons at LEP2 yield lower bounds of 
Mh,H > 92.8 GeV and Ma > 93.4 GeV. The range 0.7 < tg/3 < 2.0 in the MSSM is excluded by 
the Higgs searches for a SUSY scale Msusy = 1 TeV at the LEP2 experiments [HI. Presently 
and in the future Higgs bosons can be searched for at the Tevatron at Fermilab |6j], a proton- 
antiproton collider with a center-of-mass energy of 1.96 TeV, and the proton-proton collider 
LHC (Large Hadron Collider) with 14 TeV center-of-mass energy 01 as well as a future linear 
e^e~ collider with a center-of-mass energy up to about 1 TeV [Sj. 

The mass degeneracy of the fermions / and their superpartners / is removed by the soft 
SUSY breaking terms, which induce mixing of the current eigenstates fi and fR in addition. 
The sfermion mass matrix in the current eigenstate basis is given bj0 



- [ Ml, MlJ-\ m,{A, - y.r,) M]^ + m) ) 



with the parameters = l/r„ = tg/3 for down- and up-type sfermions. The mass eigenstates 
/i,2 of the sfermions / are related to the current eigenstates /l^/j by mixing angles 

A = fi cos 9f + fRsin 9 f 

~h = -hsmOf + fRCOsef, (2) 

which are proportional to the masses of the ordinary fermions. Thus mixing effects are only 
important for the third-generation sfermions t, b, f. The mixing angles acquire the form 



sm2ef = '\ ' ; , cos2% = ^ 1^ (3) 

Jl 72 Jl 72 



and the masses of the squark mass eigenstates are given by 

1 



2 2 I 

= TUf + 

/l,2 / 



M? + Ml T J {Ml - Mlf + Am){Af - /xr^)2j . (4) 



The topic of this paper is the calculation of the NNLO SUSY-QCD and top-induced SUSY- 
electroweak corrections to the effective bottom Yukawa couplings. These results will affect 
all processes to which the bottom Yukawa couplings contribute, i.e. in particular the Higgs 
decay widths and Higgs radiation off bottom quarks at hadron colliders which constitutes the 
dominant Higgs boson production channel for large tg/3 at the Tevatron and the LHC 



^For simplicity, the £)-terms have been absorbed in the sfermion mass parameters Af~ 

Jl/r 
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2 Effective Bottom Yukawa Couplings 



The leading parts of the SUSY-QCD and SUSY-electroweak corrections to bottom Yukawa 
couphng induced processes can be absorbed in effective bottom Yukawa couphngs. These 
contributions correspond to the hmit of heavy supersymmetric particle masses compared to the 
typical energy scales of the particular process. The accuracy of this heavy mass approximation 
has been investigated for neutral MSSM Higgs decays into bottom quarks h/H/A — )■ hb (lol |. 
charged Higgs decays to top and bottom quarks — )■ th jlli] and Higgs radiation off bottom 
quarks at e"'"e~ colliders 12] and hadron colliders 13|, [ij] by comparing it to the full NLO 
results. For large values of tg/3 the approximation turns out to agree with the NLO results to 
better than one per cent. 



2.1 Effective Lagrangian 



The leading corrections can be obtained from the effective Lagrangian [10|, lU 
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with the one-loop expressions (C^ = 4/3) [1 
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The generic function / is defined as 
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The fields 0? 02 denote the neutral components of the Higgs doublets coupling to down- 
and up-type quarks, respectively. They are related to the mass eigenstates H, A by 
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[vi + H cos a — h sin a + iA sin (3 — iG^ cos /?] 



[f 2 -|- H sin a -f- /i cos a + iA cos /3 -f sin (3] 

V2 



(8) 
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The two vacuum expectation values are related to the Fermi constant Gp, v 
1/\/2Gf- The would-be Goldstone field is absorbed by the Z boson and generates its 
longitudinal component. The top Yukawa coupling Af is related to the top mass by rrit = 
\tV2/\/2 at lowest order. The corrections A;, induce a modification of the relation between the 
bottom quark mass mb and the bottom Yukawa coupling A;,, 



mb = -j=-[l + Ab] 

The effective Lagrangian of Eq. ([5]) can be parametrized as 

^eff = -—b[ g^h + g^ H-gt tl,A]b 



(9) 



(10) 



with the effective (resummed) couplings 
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Although the SUSY corrections are loop suppressed, they turn out to be significant for 
large values of tg/9 and moderate or large values. In these cases they constitute the dominant 
supersymmetric radiative corrections to the bottom Yukawa couplings. It should be noted that 
the effective Lagrangian in Eq. ^ has been derived by integrating out the heavy SUSY particles 
so that it is not restricted to large values of tg/3 only. In order to improve the perturbative 
result it has been shown that the Lagrangian of Eq.© resums all terms of O [(a.^ tg/3)"] and 



C [(Aj tg/5)"^] The additional resummation of terms proportional to asAb [lO| will be 



neglected in this paper. 



2.2 Low Energy Theorems 

The determination of higher-order corrections to the effective bottom Yukawa couplings would 
usually require the calculation of the corresponding three-point functions. This calculation can 



be reduced to the evaluation of self-energy diagrams by the use of low energy theorems [16 
The basic idea is that any matrix element with an external Higgs boson can be related to the 
analogous matrix element without the external Higgs particle in the limit of vanishing Higgs 
momentum by the simple replacements fi — )■ a/20? ^^"^ ^2 — > v^02* ^^e latter. Thus we 
only need to compute the corresponding pieces of the bottom quark self-energy. The leading 
pieces A^ emerge from the scalar part ^^(mfc) of the self-energj|^ giving rise to the following 



^The fcrmionic self-energy can be decomposed into a scalar, vectorial and axial- vectorial part according to 
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Figure 1: One-loop diagrams of (a) the SUSY-QCD and (h) the top-induced SUSY-electroweak 
contributions to the bottom self-energy with the off-diagonal mass insertions corresponding to 
the corrections Ab of the bottom Yukawa couplings. The virtual particles involve bottom quarks 
b and squarks b, top squarks i, gluinos g and charged Higgsinos . 

relation between the pole mass nib of the bottom quark and the bottom Yukawa coupling aJI, 

rrib = + T^siruh) (12) 

where the leading terms of the self-energy ^^(mb) for large values of tg/3 are given by 

Ab = A^^^ + Af"" (13) 

The terms A^*"^ and A^'"" can be derived from off-diagonal mass insertions of the type — Ab/if 2 
in the virtual sbottom propagators and of the type \tAtV2 in the virtual stop propagators, as 
is illustrated in Figure [1] at one-loop level. The result of these diagrams is given by the finite 
expressions of Eq.([H]) after rotation of the fields in the current eigenstate basis to the mass 
eigenstates. These expressions are not renormalized since there is no tree-level coupling of 
bottom quarks to the Higgs field 02*- By means of power-counting arguments it can be proven 



that Eq. ([5]) includes and resums all leading terms of O [{as fi tg/3)"] and O [{Xf At tg/9)"] [lO|, lU 



3 NNLO Corrections 



The calculation of the NNLO corrections to the effective bottom Yukawa couplings requires 
the determination of the leading NNLO corrections to the bottom self-energy. The full NNLO 



results of the bottom quark self-energy have been presented in Refs. [17|, ll8|, ll9j before. Our 
results for the self-energy at NNLO can be extracted from their calculation by keeping only 
the terms leading in tg/3 in principle. The renormalization, however, has to be adjusted to 
our scheme in order to compare the final results. The results of Ref. 17|] are only shown for 



equal SUSY-breaking squark masses M| = M| for the full relation between the pole quark 



mass and the DR quark mass. However, since the self-energy is a symmetric function of both 



•^The bottom Yukawa coupling is not renormalized at 0{asjjXg/3) and 0{Xf At tg/S) in the effective La- 
grangian, but receives only non- leading contributions in tg/3 [10|, which do not contribute to the effective 
Lagrangian. 
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sbottom masses m^^ and m^^ their results for the linear at terms can be compared to our resuhs 
in the equal SUSY mass case. Moreover, our calculation of the bottom self-energy has been 
performed for the leading term for small bottom momenta so that we can compare directly 
with the results of Ref. |19|] which displays results in the same limit. 



3.1 Two-Loop Diagrams 

We have computed the supersymmetric QCD corrections to the effective bottom Yukawa cou- 
plings of Eq. ( ITT]) , i.e. the SUSY-QCD corrected effective Lagrangian of Eq.^S]) at NNLO. The 
two-loop calculation split up into the (9(a^//tg/9) corrections to A^'"'^^^-' and the 0{asXt^t^g(3) 
corrections to A^'"'^^''. The relevant two-loop diagrams contributing to the bottom self-energy 
are shown in Figure [21 For the determination of the NNLO corrections single scalar mass in- 
sertions in each bottom- and top-squark propagator analogous to the one-loop level have to be 
included and added in all possible ways. This means that e.g. the first diagram leads to two 
contributions, one with the mass insertion in the left sbottom propagator and the other with the 
mass insertion in the right sbottom propagator as depicted in Fig. [31 We have performed the 
whole calculation in dimensional regularization. The bottom quark momentum and its mass 
have been put to zero while keeping the bottom Yukawa coupling finite in the mass insertions 
and the charged Higgsino couplings. All supersymmetric particles as well as the top quark have 
been treated with full mass dependenc^. Since the two-loop contributions are computed in the 
limit of vanishing external momenta the calculation requires the determination of the leading 
two-loop vacuum integrals in the heavy mass expansion. These integrals can be expressed in 
terms of the one- loop one-point integral Ao(m) [20|] 

^o(™)^A-"/|^j^ (14) 



where /i denotes the 't Hooft mass, and the two-loop master integral Ti34{mi, 171^,1714) [21 



where we work in n = 4 — 2e dimensions. Integrals with higher powers of the corresponding 
propagators can be reduced to the basic integrals Aq and T134 by means of integration-by-parts 



methods [22|. In this way all singularities have been separated from the finite contributions 



analytically as poles in the parameter e. 



3.2 Renormalization 

(2) 

Contrary to the finite one-loop results, the two- loop corrections are UV- diver gent. A 
finite and /i-independent result is obtained after renormalization of the masses and couplings 

''We have neglected the tiny contributions of 0(m^{At—^J,/tgf3)/Mg^gY) of off-diagonal stop propagators 
in the gluino self-energy of the last two-loop diagram (for q = t, q = t) oi Fig. [5^ and in the gluino mass counter 
term Srrig in Eq. (1191) . 
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(a) 




Figure 2: Two-loop diagrams of (a) the SUSY-QCD and (b) the top-induced mixed SUSY- 
Q CD /electro weak contributions to the bottom self- energy involving bottom quarks h, bottom h 
and top t squarks, gluons g, gluinos g and charged Higgsinos h^. The squark-quark contributions 
to the gluino propagator have to be summed over all quark/squark flavors q/q including both 
directions of the flavor flow due to the Majorana nature of the gluino. 
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Figure 2: Cont'd. 



contributing to the one-loop result. The heavy masses mi,,mi_,mg appearing in the propagators 
have been renormalized on-shell. The trilinear coupling At has been treated in the on-shell 
scheme, too, where its counter term is derived from the counter terms of the stop mixing angle 
and masses, and the top Yukawa coupling. 
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where 6mt denotes the on-shell counter term of the top mass, 
dmt 2 f ^oim) 



rrtt 



m 



t2J 



(17) 



^Note that in the desired order the mixing angle counter term does not depend on the momentum of the 
ofF-diagonal self-energy contribution, since the gluino contribution does not contribute to our perturbative order, 



bdt = sin 20t cos 204^0(^4 J - ^o(™tJ]/(™| - m]^, 
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Figure 3: All possible mass insertions into the first diagram of Fig. contributing to the leading 
terms o/ Af, at NNLO. 

The strong coupling and the top Yukawa couphng have been defined in the MS scheme 
with 5 active fiavors, i.e. the top quark and the supersymmetric particles have been decoupled 
from the scale dependence of the strong coupling Care has to be taken to include only 

the desired order, i.e. 0{a'^,fj.tg/3 /Msusy) for A^^^ and 0{as\fAttg(3/MsusY) for Af"". In 
this order the trilinear coupling At only receives a finite renormalization at NLO. 
The complete set of counter terms is given by 




rrig + drrig 
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with the following explicit expressions {Ca = 3) 

5ml = -^Cf9s {2v4oK-J - 2Ao(m^) - 2AoK) - 4m|5o(m| ; 0, 



Oi 



Srrig 
m-n 



+2(m? - ml)Boiml ;mg, 0) 



2Eo(m?; 0, m^) - 25i(m2; 0, m^) - 



(47r) 



^ = c^^r(i + 6)(47rr^|-- + iog4| + ^c^F^7.M^iK;^5'^*J + ^iK^^ 

At 77 4 e /i J 

5q: 



TT 



^1 



:iir(l + e){A.y ^ ( -T + log ^ ) /3o + ^ log ^ + - log 3 



4^ 24 



2 

ml 



6 



mt 



A, 



mr 



-4 



ml Bo{ml ; 0, m^- J - m|Eo(m| ; 0, m 



t2' 



2 2 

mj — m| 

Cl 12 



+2- 



ti 



i2 



t2' 



2 2 

mi — mi 



tl t2 

where the lowest order beta function is given by 

„ _3Ca-Nf-1 
Po — 



(19) 



(20) 



with Np = 5 light flavors (m^ = 0) and summing over Np + 1 = 6 quark/squark flavors in the 
expressions above. The masses of all quarks except the top have been put to zero. In addition 
to the one- loop scalar integral Aq of Eq.( IT^ the residual scalar one- loop integrals are defined 



as 



20 
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(21) 
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A few comments are in order. Since we calculate the 0{a'^fitg(3) and 0{asXtAttgf3) corrections 
we only included contributions to these orders in our counter terms. This plays a role for the 
counter terms of the gluino mass, the stop masses, the top Yukawa coupling and the trilinear 
coupling At. In the gluino mass counter term we have neglected contributions of the order 
mt/rrig, in the counter terms of the stop masses and the top Yukawa coupling Xf we deleted 
terms of 0{mtAt), and finally in the counter term of the trilinear coupling At we only included 
contributions linear in At. Moreover, we have put the bottom mass rub to zero everywhere 
apart from the Yukawa coupling Xf, in the mass insertions and the Higgsino couplings. This 
leads to the simplified counter terms listed above [2^ . 

The counter term for the strong coupling constant contains the logarithmic finite terms 
required for the decoupling of the top quark and the supersymmetric particles from its scale 
dependence, i.e. the coupling as runs with only five light flavors. 



with (Np = 5) 



11Ca-2Nf 
12 



(22) 



(23) 



Anomalous Counter Terms. A complication occurs in supersymmetric theories if dimensional 
regularization is employed, since in n 7^ 4 dimensions the gluons acquire n—2 degrees of freedom, 
while the gluinos still possess 2 degrees of freedom. Thus a mismatch between the gluons and 
their supersymmetric partners is introduced at 0{e) which violates the supersymmetric relations 
between related couplings and masses. In order to restore supersymmetry at the renormalized 
level finite anomalous counter terms have to be introduced which, however, are uniquely fixed 



by the corresponding supersymmetric identities [24 



In our calculation this mismatch of degrees of freedom between supersymmetric partners 
induces finite differences between the Yukawa coupling Qs of quarks, squarks and gluinos and 
the corresponding gauge coupling gg, 



9s = 9s 



1 _|_ _ ^ 
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(24) 



and between the quark Yukawa couplings (generically denoted by Xnqq), the corresponding 
Higgs couplings to squarks (generically Xngq) and the Higgsino Yukawa couplings to quarks 
and squarks (generically Xfj^-), 
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(25) 



In our calculation we have expressed all couplings in terms of the Standard Model MS Yukawa 
coupling Xnqq- The anomalous counter terms for the Hqq vertex has already been included 
in the counter term 6At given above. The left-over anomalous counter terms are those for 
the strong coupling gs arising in the SUSY-QCD corrections A?^^^^^ at the one-loop level. 
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the bottom Yukawa coupling of the mass insertions of A^'"'^^^'', which corresponds to the 
generic couphng Xnqg, the bottom Yukawa couphng A;, corresponding to A^^- generically in the 

SUSY-electroweak corrections A^'"'^^'' and the top Yukawa couphng factor Xf of A^^"'^^^ which 
corresponds to the generic product XHqqXfiqq- This resuhs in the following total anomalous 



counter terms 



23 



b;anom Ig 4 A TX ^ 13 2 I Tl 



5^tno^ = ( - - ^C. ) ^A^^^) = -C.^A^^^) (26) 



These counter terms have been added to the final results after including the counter terms of 
Eq.f lT9|) to obtain consistent final results. 

3.3 Resummation 



According to the power-counting arguments of Ref. [10| the leading NLO contributions of 
0{^tgf3) can only emerge from single mass insertions as depicted in Fig. [H i.e. NNLO terms of 
0{fiHg'^f3) are absent, since a second mass insertion as required by the second power in /itg/3 is 
suppressed by an additional power of rrih/MsusY- In complete analogy the absence of C(A^tg^/3) 
in the SUSY-electroweak part A^'"' can be proven. Thus the leading two-loop corrections are 
of (9(a^/itg/3) and 0{asXtAttg(3) respectively, arising from single mass insertions in the two- 
loop bottom self-energy diagrams contributing to T,s{mb) in Eq. f|T2l) . These are indeed all 
contributions taken into account in this work. At higher orders terms of (9(as/i^tg^/3) and 
0{asA'^tg'^(3) are absent so that our results of 0(a^/itg/3) and O^agX'fAttgP) are automatically 
resummed if included in A^ according to Eq. ([5]) with 

^QCD ^ ^QCD{1) _^ ^QCD{2) 

^dw ^ ^elwil) ^ ^elu,(2) ^27) 

where A^'"'^ '■^^ and A^'™ denote our novel NNLO corrections. These expressions thus resum 
all terms of C(a>"tg"/3), C(a^+V"tg"/3), C(A2"AJ'tg"/3) and C(a,A2"Aj^tg"/3). 

3.4 Limits 

The final results for the two-loop corrections A^*"^ and A^'"' "-^^ are too lengthy to be displayed 
here. However, they can be given explicitly in certain limits: 
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A eiu;(l) 



TT 

i2 



2 (47r)2 M 



- ^^vl4 + 2^°^;;a^/^^ ^^^^ 

where terms of 0{ml / M"^) have been neglected. The logarithmic top mass term of A^*"'^*-^'* 
can be absorbed in the running strong coupling by evolving it with 6 active flavors. However, 
since the top mass is not much lighter than the supersymmetric masses this logarithm does not 
become that large that this change of scheme is required. The last logarithm compensates the 
scale dependence of the strong coupling ^^(/i/j) in the one- loop expression We have 

verified that our results for A^*"^ can be obtained from the Xi, terms of Eq. (4.8) in Ref. 19 



if the scheme transformations and decoupling relations of the strong coupling and the trilinear 
coupling Ah are taken into account properly. Moreover, we have found agreement with the 



terms of Eq. (61) in the first paper of Ref. 17(1 after appropriate scheme transformations of 
as, rUg, mi^^^ and the bottom mass and the necessary subtraction of the product of the vectorial 
part of the one-loop self-energy Ey(mb) and the scalar part T,s{mi,) which contributes to the 
relation between the pole quark mass and the DR mass in addition to the pure self-energy at 
the two-loop level. 

The final expressions for A^''" exhibits a logarithmic term of the common supersymmetric 
mass M which can be absorbed in the low-energy Yukawa couplings of the charged Higgsinos 
to stops and bottom quarks which have to be defined with properly decoupled supersymmetric 
contributions. These low-energy Yukawa couplings have to be evaluated at scales of the order 
of the top quark mass. This, however, is only necessary if the mass splitting between the top 
quark and the supersymmetric particles becomes large, giving rise to logarithms of their mass 
ratios. This does not occur in realistic scenarios so that we did not perform this low-energy 
decoupling. 



^It should be noted that degenerate squark masses require A}, — /itg/3 and At = /i/tg/3. Moreover, due to 
the different Z?-terni contributions to the 5[/(2)-relation between the left-handed stop and sbottom mass terms 
the equality of the stop and sbottom masses cannot be realized. However, we display the full result in this limit 
as a good approximation for small mass splittings. 
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The next non- leading term in this limit is of 0(log ^ m?/m? ) for the SUSY-QCD part A?'"'^ 

and of (9 (max{ /x^, m|}/m?) for the SUSY-electroweak part A^'"'^^''. The expression of A^'"'''^'* 
is not altered from Eq. ([6]) in this limit so that we do not display it again. The logarithms of the 
top and squark masses in A^'"'^ ^'^^ can be absorbed in the running QCD coupling as if the latter 
includes the top quark and all squarks in its /3 function. Since the real mass splittings are not 
so large that these logarithms become sizeable this procedure is not necessary. Again we have 
verified that our results for A^*"^ in this limit can be derived from the Xf, terms of Eq. (4.10) 
in Ref. [l^ by taking into account the decoupling properties and scheme transformations of as 
and Ab consistently and identifying all squark masses. As for the previous limit the logarithm of 
the gluino mass in A^'"" can be absorbed in the low-energy Yukawa couplings of the Higgsinos 
to stops and bottom quarks by decoupling the gluino contributions consistently. Since these 
logarithms do not become large in realistic cases we did not perform this decoupling but kept 
the full supersymmetric relations to the related Yukawa and gauge couplings. 



4 Results 

The results obtained in this work have been inserted in the program HDECAY [25|] which 



calculates the MSSM Higgs masses and couphngs according to Ref. [26| as well as all partial 
decay widths and branching ratios including the relevant higher-order corrections . For large 
values of tg/3 the decays of the neutral Higgs bosons are dominated by the decays into bb and 
r+r". Their branchin g ra tios had been studied with the one- loop expressions of the correction 



A;, of Eq. in Ref. [10[- It has been demonstrated that the scale dependence of the strong 
coupling as{f^R) of Eq. ([6]) implies a theoretical uncertainty up to the 10% level for scenarios 
of large SUSY-QCD corrections to the bottom Yukawa couplings such as, e.g., the 'small cte//' 
scenario [2j 
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4.1 Higgs Decays into Bottom Quark Pairs 



The partial decay widths of the neutral Higgs bosons $ = h,H,A into bottom quark pairs, 
including QCD and SUSY-QCD corrections, can be cast into the form l3 



r[$ ^ bb] 



4^2 



-m, 



n 



:(M^) [Aqcd + Af ] gf [gf + A— ^] 



(30) 



where regular quark mass effects are neglected. The large logarithmic part of the QCD cor- 
rections has been absorbed in the running MS bottom quark mass mf,(M<j,) at the scale of 
the corresponding Higgs mass M$. The QCD corrections Aqcd and the top quark induced 



contributions Af read as 



28 



A 



QCD 



1 + 5.67^^ + (35.94 - 1.36iV,) 



TT 



IT 



+ (164.14 - 25.77Nf + 0.259Nj.) 

2 



TT 



A 



h/H 



h/H , 
9t I 



h/H 



9b 

gf fas{MA) 



(MvH) y 



TT 



rnl{Mh/H) 



h/H 



TT 



Ml 



mi{MA 
Ml 



(31) 



where Np = 5 active flavors are taken into account. In the intermediate and large Higgs 
mass regimes the QCD corrections reduce the bb decay widths by about 50% due to the large 
logarithmic contributions. 



The main part of the SUSY-QCD corrections [29| has been absorbed in the effective bottom 
Yukawa couplings gf of Eq. (fTTj) . The remainder A^'g^^ is small in phenomenologically relevant 
scenarios for large values of tg/9 10|. It should be noted that this remains true even for very 
light SUSY masses in the range of 0(100 GeV). This can be understood by deriving the 
asymptotic expression of A^Fg^^ for large values of the Higgs masses compared to the SUSY 
masses (we neglect the bottom mass here). 



A rem , ^ CX-sj^-Il) 

^SQCD '-'F — — mg 



TT 



(A - r/^/itg/3) I{ml^,ml^,ml) 



(32) 



with the coefficients 



Vh 



tga 



1 



VA 



tg/3 tg2/3 

for the heavy scalar H and the pseudoscalar A. Thus the large Higgs mass limit approaches a 
finite expression of 0{Ab/ iitg(3) relative to the leading A;, terms for large values of tg/3. This 
explains the validity of the Af, approximation even for small SUSY masses compared to the 
Higgs masses. In the following we include the mixed top Yukawa coupling induced SUSY- 
QCD/electroweak corrections in the couplings gf, too. 
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4.2 Numerical Results 



The numerical analysis of the neutral Higgs boson decays into bottom quark pairs is performed 
for two MSSM benchmark scenarios 27] as representative case^: 



small cteff- 



tg/3 

M2 : 



: 30, Mq = 800 GeV, A'% = 1000 GeV, 
500 GeV, Ah = At = -1.133 TeV, /i = 2 TeV 



gluophobic: tg/3 = 30, Mq = 350 GeV, Mg = 500 GeV, 

M2 = 300 GeV, Ab = At = -760 GeV, = 300 GeV (33) 

We use the RG- improved two- loop expressions for the Higgs masses and couplings of Ref . jiij . 
Thus the leading one- and two-loop corrections have been included in the effective mixing angle 
a. The top pole mass has been taken as = 172.6 GeV, while the bottom quark pole mass 
has been chosen to be rrib = 4.60 GeV, which corresponds to a MS mass rrihi^b) = 4.26 GeV. 
The strong coupling constant has been normalized to as{Mz) = 0.118. We have used the top 
and bottom pole masses in the mass matrices of the stop and sbottom states as implemented 
in HDECAY. 

In the following we will not address the parametric uncertainties of the corrections to the 
bottom Yukawa couplings and branching ratios, which will be significantly reduced after the 
LHC and a linear e'^e~ collider with energies in the range up to about 1 TeV [8]. However, we 
will concentrate mainly on the renormalization scale dependence of our final results as a first 
estimate of the residual theoretical uncertainties after including our novel NNLO corrections 
to the bottom Yukawa coupling^. In Fig. H] the scale dependence of the A^*"^ and A^'*" 
terms is shown for the small Q!e// scenario and in Fig. E] for the gluophobic scenario. The 
^QCD (.QTj.Ti.Q(.i\QT^ amounts to about 90% in the small Oe// scenario partially compensated by 
a negative /S.f'^ contribution of C(10%). By comparing the dashed one-loop bands with the 
full two-loop curves a significant reduction of the residual scale dependence can be observed 
leading to a final scale uncertainty in the per-cent range at NNLO. An analogous reduction 
of the scale dependence can be observed in the gluophobic scenario, where the /s!^^^ terms 
amount to about 35% and are significantly reduced by about 20% due to the electroweak Af^ 
contribution. In general, however, the sizes and signs of these contributions strongly depend on 
the MSSM scenario. Moreover, broad maxima/minima develop at about 1/3 of the central scales 
determined by the average of the corresponding SUSY- masses, i.e. /iq = (m^^ + m^^ -|- mg)/3 
for A^*""^ and fiQ = {mi^ + m^^ + /x)/3 for /S.f'^ at NNLO in contrast to the monotonic scale 
dependence of the one-loop expressions. For the central scale choices the two-loop corrections 

^Note that in the small cte// scenario we increased the gluino mass by a factor of two relative to Ref. [l^] in 
order to enhance the size of Af,. This value for the gluino mass has also been used in Ref. [2^ so that Eq. (5) 
of that work has to be corrected. A general analysis of the impact of the Af, terms in the MSSM is beyond the 
scope of our work so that we restrict our numerical analysis to these two scenarios. 

^The scheme and scale choices of the stop, sbottom and gluino masses as well as the trilinear coupling At 
modify the results by a few per cent at NLO already so that their impact ranges below the per-cent level at 
NNLO. These effects are neglected in our analysis, since the dominant uncertainties originate from the scale 
choices of the strong coupling constant and the top Yukawa coupling At. 
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to A^*^'^ amount up to (9(10%), while they are in the per-cent range for Af"' 0. Since the 
corrections A;, enter the effective Yukawa couphngs as in Eq. f|TT]) they have an immediate 
impact on all processes induced by these bottom Yukawa couplings. 

As a particular example we show the partial decay widths of the neutral MSSM Higgs 
bosons into bb pairs in Figs. E] and [7] in the small Oe// and the gluophobic scenario, respectively. 
The two-loop corrections to A;, lead to negative corrections to the partial decay widths for 
the central scale choices of C(10%). The bands at one-loop order (dashed blue curves) and 
two-loop order (full red curves) are defined by varying the renormalization scales of A^'"'^ and 
Af^ independently between 1/2 and 2 times the corresponding central scales Hq. A significant 
reduction of the dashed one-loop bands of C(10%) to the full two-loop bands at the per-cent 
level can be inferred from these results. The small kinks in the partial decay width r{H bb) 
in the gluophobic scenario originate from the bibi and the 62^2 thresholds in A^'g^^ of Eq. (130|) . 
The small kink in the partial decay width r(y4 — )■ bb) on the other hand can be traced back to 
the 6162 and 62^1 thresholds. The light scalar Higgs boson decay width r{h — )■ bb) exhibits a 
strong dependence on the pseudoscalar mass Ma for Ma ^ 120 GeV, since it slowly approaches 
the SM-limit for large Ma- The drop of r{h — )■ bb) for Ma ~ 145 GeV in Fig. [H]is due to the 
vanishing of the effective mixing angle a in the MSSM Higgs sector. 

The scale uncertainties in the partial decay widths translate into theoretical uncertainties in 
the corresponding branching ratios of the neutral MSSM Higgs bosons. These errors cancel to 
a large extent in the dominant branching ratio BR{^ — )■ bb), but induce sizable uncertainties of 
the non-leading branching ratios at one-loop order. The branching ratios for the two dominant 
decay modes into bb and r+r^ pairs are depicted in Figs. [S] and M for the small aeff and 
the gluophobic scenario, respectively. The uncertainties of the branching ratios reduce from 
O (10%) at one- loop order to the per-cent level at two-loop order. The per-cent accuracy now 
matches the expected experimental accuracies at a future linear e^e~ collider. The thresholds 
visible in the branching ratios of the heavy scalar Higgs boson in Fig. IHlare due to the kinematical 
opening of the titi, bibi and 62^2 decay modes in consecutive order with rising Higgs mass. The 
thresholds of the pseudoscalar branching ratios are due to the decay modes into X1X25 XiXi ^iih 
and 6162- 

Since we have determined the effective resummed bottom Yukawa couplings at two-loop 
order the results will also affect all other processes which are significantly induced by bottom 



Yukawa couplings, as e.g . MSSM Higgs radiation off bottom quarks at e~^e colliders [12|, |30 



and hadron colliders j3l| . The two- loop corrections can easily be included in the corresponding 



numerical programs. However, these analyses are beyond the scope of this work. 



5 Conclusions 

In this work we have determined the NNLO corrections to the effective bottom quark Yukawa 
couplings within the MSSM for large values of tg/3, since the theoretical uncertainties are 
sizable at NLO in these regimes. The leading parts of the SUSY-QCD and top-induced SUSY- 
electroweak corrections originate from factorizable contributions due to virtual squark, gluino 
and Higgsino exchange, which can be absorbed in effective Yukawa couplings in a universal 
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Figure 4: Scale dependence of the SUSY-QCD correction a'^^'^ and the SUSY-electroweak 
correction Af'^ at one-loop and two-loop order in the small a^ff scenario. 
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Figure 5: Scale dependence of the SUSY-QCD correction A^'"'^ and the SUSY-electroweak 
correction Af'^ at one-loop and two-loop order in the gluophobic scenario. 
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Figure 6: Partial decay widths of the light scalar h, the heavy scalar H and the pseudoscalar A 
Higgs bosons to bb in the small aeff scenario. The dashed blue bands indicate the scale depen- 
dence at one-loop order and the full red bands at two-loop order by varying the renormalization 
scales of A'^'"^ and Af'^ independently between 1/2 and 2 times the central scale given by the 
corresponding average of the SUSY-particle masses. 
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Figure 7: Partial decay widths of the light scalar h, the heavy scalar H and the pseudoscalar A 
Higgs bosons to bb in the gluophobic scenario. The dashed blue bands indicate the scale depen- 
dence at one-loop order and the full red bands at two-loop order by varying the renormalization 
scales of A'^'"^ and Af"^ independently between 1/2 and 2 times the central scale given by the 
corresponding average of the SUSY-particle masses. 
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Higgs bosons to bb and t^t~ in the small Oe// scenario. The dashed blue bands indicate the 
scale dependence at one-loop order and the full red bands at two-loop order by varying the 
renormalization scales of A'^'"^ and Af"^ independently between 1/2 and 2 times the central 
scale given by the corresponding average of the SUSY-particle masses. 
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way. We have calculated the two-loop SUSY-QCD corrections to these effective bottom Yukawa 
couplings. 

In summary, the significant scale dependence of O (10%) of the NLO predictions for pro- 
cesses involving the bottom quark Yukawa couplings of supersymmetric Higgs bosons require 
the inclusion of NNLO corrections. For the corrected Yukawa couplings, we find a reduction of 
the scale dependence to the per-cent level at the NNLO level. The improved NNLO predictions 
for the bottom Yukawa couplings can thus be taken as a base for experimental analyses at the 
Tevatron and the LHC as well as the ILC. 



Note added 

During the publication procedure of our work another paper appeared 32| where the bottom 
Yukawa couplings have been calculated at two-loop order within SUSY-QCD. After adjusting 
the different renormalization schemes and approximations we found mutual agreement between 
the two calculations. 
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